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We present muon-spin rotation measurements on polycrystalline samples of the complete family of 
the antiferromagnetic (AF) zigzag chain compounds, Na^Cai-a; V2O4. In this family, we explore the 
magnetic properties from the metallic NaV2 04 to the insulating CaV^CU. We find a critical x c (~ 
0.833) which separates the low and high Na-concentration dependent transition temperature and its 
magnetic ground state. In the x < x c compounds, the magnetic ordered phase is characterized by a 
single homogenous phase and the formation of incommensurate spin-density-wave order. Whereas 
in the x > x c compounds, multiple sub-phases appear with temperature and x. Based on the muon 
data obtained in zero external magnetic field, a careful dipolar field simulation was able to reproduce 
the muon behavior and indicates a modulated helical incommensurate spin structure of the metallic 
AF phase. The incommensurate modulation period obtained by the simulation agrees with that 
determined by neutron diffraction. 

PACS numbers: 76.75.-K, 75.50.Lk, 75.50.Ee 



I. INTRODUCTION 



Experimental and theoretical studies of quasi-one di- 
mensional (quasi-lD) magnets have flourished recently. 
On theoretical grounds, such systems seem to experi- 
ence low dimensionality combined with frustrated inter- 
actions of nearest-neighbors and next-nearest-neighbors 
competition, leading to spin density waves forming chi- 
ral order [l| and other exotic ground states [2, 3]. Ex- 
perimentally, progress has accelerated due to recent use 
of high-pressure techniques UHl], enabling synthesis of 
numerous compounds. Such materials show remarkable 
physical behavior @-|||, where the nature behind these 
phenomena is governed by a strong spin-spin interaction 
along one direction combined with a much weaker inter- 
action in the other directions. A novel quasi-lD spin sys- 
tem of particular interest is NaV2 O4 (NVO ) (lpj . which 
is isostructural with CaV2 04 (CVO)[Tl[. Unlike other 
vanadate compounds, which are spinels [l2j. NVO and 
CVO posses an orthorhombic structure with the Prima 
space group, where irregular hexagonal ID channels are 
formed by a series of edge-sharing VOg octahedra aligned 
along the 6-axis. The magnetic V ions are thus ar- 
ranged in a ID zigzag chain. Despite their closely related 
chemical and crystallographical nature, they exhibit dra- 
matically different magnetic and electronic properties. 
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The metallic NVO incorporates a mixed valence V +3 5 
state and undergoes an antiferromagnetic (AF) transi- 
tion at Tn = 140 K, with a complex magnetic ground 
state[l(J[l3j]. In fact, recent neutron studies indicate the 
formation of incommensurate magnetic order in NVO be- 
low T]v[lJ]. The insulating CVO, on the other hand, 
has a V +3 valence {t\ g ,S = 1) with an AF transition 
at Tn = 70 K, coexisting with intrachain ferromagnetic 
interactions. Furthermore, the chemical substitution of 
Na for Ca causes an interesting insulator to metal tran- 
sition at a critical Na concentration x > x c ~ 0.83 in 
Na x Cai_ 2 ;V204 (NCVO). The presence of the insulator 
to metal transition suggests a dynamic change in the elec- 
tronic structure of NCVO with x, resulting in the change 
in its magnetic ground state. It is therefore desirable 
to systematically investigate the evolution of the micro- 
scopic magnetic order/disorder in NCVO. 

Recently, through the use of positive muon-spin rota- 
tion and relaxation (/i + SR), we found the formation of 
static AF order for a few NCVO compounds below their 
Tn and proposed a magnetic phase diagram[l3|]. /i+SR 
is known to be a powerful technique for studies on mi- 
croscopic magnetic nature in solids, due to its inherent 
100% spin polarization and local sensitivity. The most 
interesting feature in NCVO is the unusual coexistence 
of AF order and metallic behaviour below Tn for NCVO 
with x > x cl: . The main purpose of this paper is, there- 
fore, to clarify the microscopic magnetic ground state 
of the entire NCVO family through the use of /i + SR and 
numerical simulations based on dipolar field calculations. 
This is, to our knowledge, a pioneering attempt to deduce 
the incommensurate modulation period from /i + SR data. 
We also seek to complete the magnetic phase diagram of 
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NCVO by additional /i + SR measurements on insulating 
and metallic NCVO samples with x = 0.166, 0.41, 0.875 
and 0.958, which were not described in our previous re- 
port. Our results reveal that the magnetic ground state 
in the insulating compounds with < x < also 
an incommensurate magnetically ordered phase. This 
means that not only pure CaV204 but also the complete 
family of NCVO exhibit incommensurate magnetic order 
at low T, concurrent with the change of the valence state 
of the V ion from 3 to 3.5 with increasing x. This could 
provide a very interesting challenge for theorists: to pre- 
dict the magnetic phase diagram of this zigzag chain 
compounds. 
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II. EXPERIMENT 

Polycrystalline samples of NCVO were prepared by a 
solid-state reaction technique under a pressure of 6 GPa 
using CaV204, Na4V20y, and V2O3 powders as starting 
materials. A mixture of the three powders was packed in 
an Au capsule, then heated at 1300°C for 1 hour, and fi- 
nally quenched to ambient T. A powder X-ray diffraction 
(XRD) analysis showed that all the samples were almost 
single phase with an orthorhombic system, Pnma space 
group, at ambient T. dc-% measurements showed that all 
our samples have almost the same T dependence as that 
in the previous report [l5[ . The preparation and charac- 
terization of the samples have been reported in greater 
detail elsewhere [l5| . 

The ^ + SR measurements were carried out on the M20 
surface muon beamline with the LAMPF spectrometer 
at TRIUMF, the Canada's National Laboratory for Par- 
ticle and Nuclear Physics located in Vancouver, Canada, 
which provides a highly intense, polarized beam of pos- 
itive muons. The samples, 13 polycrystalline discs with 
6 mm diameter and 5 mm thickness of NCVO, were 
mounted on the sample holder for the muon veto cryostat 
insert [16| by Mylar-tape. 

Weak transverse field /i + SR (wTF-/i + SR) and zero 
field /z + SR (ZF-^ + SR) spectra were measured in a non- 
spin-rotated mode- i.e., the initial muon-spin direction 
is anti-parallel to the muon momentum direction. Here, 
"weak" means that the applied field is significantly less 
than any possible spontaneous internal fields (H- m t) in 
the ordered state. A wTF-/z + SR technique is sensitive 
to local magnetic order through the /i + spin polarization 
amplitude. In contrast, ZF-^ + SR is a sensitive probe of 
local magnetic [dis] order through the precession of the 
muon due to internal magnetic fields at the muon inter- 
stitial sites. 



III. RESULTS 

In the inset of Fig. [T^b) we show the typical wTF 
asymmetry (^4tf) represented here by the x = 0.41 sam- 
ple at two temperatures, above and below Tn. In the 
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FIG. 1: (Color online) (a) Temperature dependence of the 
normalized TF asymmetry, iV^TF, for the complete x family, 
(b) The transition temperature, Tn versus the Na concentra- 
tion, x. Inset, a typical wTF-^t + SR spectrum of the NCVO 
samples. Shown here the paramagnetic phase, T — 85, and 
the ordered state T — 1.8 K of the x = 0.41 sample. The 
solid lines indicates a fit to Eq. (jTJ) . 



paramagnetic state, at high temperatures (T > 80 K), 
the muon-spin precess at the frequency corresponding to 
the external transverse field Htf = 50 Oe. Below Tn, 
Atf decreases, due to the strong static internal mag- 
netic fields, Hi n t, obeying H m t > -Htf- The remaining 
oscillating amplitude, shown in the figure, reflects the 
portion of muons not coupled to -ffmt, thus the fraction 
of the wTF oscillatory signal below Tn is negligibly small 
for the magnetic nature of these compounds. The wTF 
spectrum is well described by an exponentially relaxing 
cosine oscillation; 



A P TF (t) = A TF exp[-(A T F<0] cos(uj TF t + 0), (1) 

where wtf = 7^-ffTF with 7 /1 /(27r) = 13.554 kHz/Oc. 
The fit is represented by the solid lines in the inset of 
Fig. lib). 

FigureQJa) depicts the temperature dependence of the 
normalized ^4tf, NAtf, which is proportional to the 
fraction of the paramagnetic regions in the sample. In 
NCVO, we find two distinct AF phases; one occurs as 
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FIG. 2: (Color online) The apodized ZF-^+SR FT of the 
x = 0.41 sample. Inset shows the time-spectrum the solid 
line represents a fit to Eq. ([3j . 



Na concentration increases up to x ~ 0.8 and the transi- 
tion temperature, Tn decreases, the second occurs with 
x > 0.833, as the Na concentration increases Tn takes a 
sharp turn and increases to higher and higher tempera- 
tures with increasing x above 100 K. These observations 
are summarized in Fig.[IJb) which plots T/v versus the Na 
concentration x. T/v is found at the middle of the transi- 
tion and the errors on T/v are taken as 10% of the tran- 
sition temperature width. Since the 0.083 < x < 0.33 
and x — 0.775 samples do not experience a sharp mag- 
netic transition, we find these samples to be a mixture 
of phases. Additionally, we can indicate that a critical 
doping level x CT ~ 0.83 separates the two phases. How- 
ever, we wish to emphasize that both bulk magnetiza- 
tion measurements and x-ray diffraction analysis show 
that each of these samples is almost a single phase with 
an orthorhombic structure [il fl3|. This clearly indicates 
the importance of the current microscopic measurements, 
which show results quite different from the macroscopic 
bulk measurements. 



A. insulating region; NCVO with < x < 0.775 

We now explore the two AF phases, identified by 
the wTF-/z + SR measurements, using the ZF-/i + SR tech- 
nique. The ZF-/x + SR is a sensitive site-based probe of 
static magnetism, in which the time dependent muon po- 
larization signal is determined exclusively by the iTnt's 
in the sample. We expect that different magnetic phases 
would cause different H- lnt s, resulting in different muon 
frequencies. Therefore, following upon the wTF mea- 
surements, we performed detailed ZF measurements for 
the entire x range. In the 0.41 < x < 0.66 range, we 
concentrate on the low T behaviour, where the AF tran- 
sition occurs for T < 10 K. In contrast, for the x > 0.833 
samples, we examine a broad T range below 150 K. 

Before discussing our current results in NCVO, we 
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FIG. 3: (Color online) (a) The ZF SDW frequency [/sdw = 
wsDw/(27r)] versus temperature and (b) the normalized /sdw, 
/iv, versus the reduced temperature [see text]. The solid line 
indicates the T dependence of the BCS gap energy. 

briefly summarize the ZF-^( + SR results of pure CVO[l3[. 
The CVO ZF-/x + SR spectra indicated four temperature 
dependent frequencies below the AF transition, two of 
which have significant amplitudes. In order to clarify the 
ground state, dipolar field calculations suggested that AF 
order exists along the two ID legs in the zigzag chain, 
and AF order between the zigzag chains. This result was 
found to be consistent with neutron measurements both 
for iT n t's and the magnitude of the ordered V moment. 

The typical raw ZF-spectrum, represented by the x — 
0.41 sample, taken at T = 1.8 K is depict in the inset of 
Fig. [2j One can clearly see a strongly damped oscillation 
in the early time domain, t < 0.15 /zs, at the lowest 
T measured. In order to fit the spectrum, we, at first, 
attempted to apply a combination of an exponentially 
relaxing cosine signal and a slowly relaxing signal. The 
former corresponds to the static AF ordered signal, and 
the latter does to the "1/3" tail for the powder sample: 

A P ZF (t) = A AF e^p(~X AF t)cos(uJAFt + (t)AF) 

+ A tail exp(-A t aiit)- (2) 

Although the ZF-spectrum was reasonably fitted by 
Eq. @, the initial phase (0af) delays by (66 ± 2)°, 
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(57 ± 9)° and (44 ± 6)° in the x = 0.41, 0.5 and 0.66 
samples. Since 0af should be zero for a simple commen- 
surate AF order, the delay indicates either a wide dis- 
tribution of Hi n t or a formation of an incommensurate 
(IC) AF order. In fact, the Fourier Transform (FT) of 
the ZF time-spectrum in 0.41 < x < 0.66 clearly demon- 
strates the wide field distribution at the muon sites [see 
the main panel of Fig. [2] as an example of the x — 0.41 
sample] . The FT spectrum also indicates the presence of 
a shoulder on the higher frequency side, suggesting the 
formation of complex magnetic order. Indeed, the ZF- 
spectrum is better fitted by a zeroth-order Bessel func- 
tion of the first kind [Jo(w M t)], especially in explaining the 
fast relaxing behaviour in the early time domain (before 
0.02 /is). A ZF-spectrum described by Jo is well estab- 
lished signature that the /x + 's experience an IC magnetic 
field in the lattice [17[. Hence, the ZF-^ + SR spectra at 
low T were fitted with a combination of two signals; 

A P ZF (t) = A S Dwexp(-A S Dwi)Jo(^SDwO 

+ A ta iiexp(-A ta ut). (3) 

The fit is demonstrated by the solid line in the inset 
of Fig. [5] In Fig. EJa), we plot the frequency /sdw = 
wsdw/27t, extracted from the fits, versus T. For all 
x < 0.66, as T increases, /sdw decreases with an in- 
creasing slope (d/sDw/dT) and approaches zero at T^. 
Interestingly, as indicated in Fig. 3b, the normalized 
/sdw [/n = /sdw (T) //sdw (T ->• K)] and the reduced 
T (T r = T/T N ) are independent of x [see FigEfb)]. In 
fact, the wn-vs.-TJ. curves now collapse to a single univer- 
sal curve, indicating the same origin of the IC transition. 
This universal curve is well explained by the T depen- 
dence of the BCS gap energy, as expected for the order 
parameter of the IC-AF state pjfl9|. This was recently 
found by a 51 V NMR study, which suggested an IC AFM 
for x > 0.77[2(| (Note that AF resonance lines in the 
metallic phase were observed in Ref.[20]). We therefore 
find that the IC phase in the insulating region includes 
all compounds, at low T. 

The transverse relaxation rate, AsdWi shows typical 
critical behavior as T — > Tjv, namely, Asdw increases 
with increasing T up to Tjv as expected [2ll423| . The 
longitudinal relaxation rate, Ataii, can in principle reveal 
whether frustration prevents the system from reaching 
the full static limit, as observed in some highly frustrated 
systems: normally one expects Ataii — ► (static limit) as 
T — >• 0. However, Ataii is very small (< 1/is -1 ) compared 
with Asdw and difficult to determine reliably, especially 
in the presence of several signals, hence we cannot make 
a confident statement. 



B. metallic region; NCVO with x > 0.83 

The FT of the whole x range, at T = 2 K, is pre- 
sented in Fig. 01 Clearly, the ZF measurements of the 
metallic compounds, i.e., x > 0.83 samples, reveal differ- 
ent behavior than those at smaller x. Unlike the insulat- 
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FIG. 4: (Color online)The apodized ZF FT of the whole x 
range, top displays the insulating region characterized by a 
single wide frequency distribution, bottom shows the metallic 
region characterized by numerous frequencies. 



ing compounds, the metallic compounds, display multiple 
frequencies (two in x — 0.875, four in x > 0.9), forming 
one frequency at higher T and vanishing above the phase 
transition (not shown). Therefore, the ZF-spectra for 
these compounds were fitted by a sum of four oscillating 
signals with a slow exponential relaxation, 



A P Z F(t) = A t cos(uit + <f>j) exp(-Xjt) 



i=l 



Ataii exp(-Ataiii) 



(4) 



where n = 4 and uii,cj>i, and Xi are the angular frequency, 
the initial phase, and the relaxation rate of the i-th muon 
precession, Ataii an d Ataii are the asymmetry and the 
relaxation rate for the "1/3 tail" signal [241 [25[, as we 
measured powder samples. 

The quality of the fit, shown for NVO, is represented 
by the solid green line in Fig. [Ha). The T dependence 
of the frequencies, fi = Wi/ (2tt), of the 0.83 < x < 0.95 
samples shown in Fig. [HJa). In the x = 0.83 sample, we 
note that at T < 70 K two frequencies are distinguished, 
but at 70 < T < 110 K a single frequency is observed. 
In x — 0.91 [x — 0.95), the number of frequencies is de- 
creased from 4, at T < 70K to 2 (3) and then to one 
frequency. Recent neutron scattering measurements sug- 
gests that a small displacement of 2 ~ ions, which leads 
to a change in the muon sites [l4[ is the trigger to this 
intermediate phase. This scenario was also proposed in 
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FIG. 5: (Color online) (a) The temperature dependence of 
the frequencies in the metallic phase, top to bottom, x = 
0.875, 0.91, 0.95. (b) The the reduced-temperature (T r ) de- 
pendence of the normalized /2 frequency, /n- The solid line 
demonstrates the BCS gap enegy. 



other compounds [26j. One of the distinctive features in 
this x range, are two common frequencies. At base T, 
/2 ~ 35MHz and f± ~ 17MHz are x independent. As in 
the insulating phase, by plotting the frequency versus the 
reduced temperature results in the collapse of these fre- 
quencies to a single curve, suggesting the IC AFM in this 
phase as well. Fig. [5^b) demonstrates this for f%. In the 
next section we aim to find the correct spin configuration 
which produces this behavior. 



IV. DISCUSSION 
A. Incommensurate AF structure for NVO 

We now address the origin of the 4 frequencies ob- 
served in metallic phase using computer simulations. We 
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FIG. 6: (Color online) (a) The raw ZF-/X+SR data of NaV 2 4 . 
The solid (grey) line indicate a fit, the helical simulation is 
represented by the dotted line (red), and the linear by the 
dashed (green) line, (b) The suggested magnetic structure, 
the moment orientation in the dc plane (right), the moment 
size and the helical modulation length I (center), the magnetic 
V ions and their moments in I = 1. 
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k[b] (A- 1 ) 


x 2 


Helical 


0.916 


0.381 


1.591 




0.958 


0.363 


3.732 




1 


0.191 


0.910 


Linear 


1 


0.571 


105.203 



TABLE I: The helical simulation result showing the modu- 



lation length along 6, 
concentrations x. 



I, and the goodness-of-fit, x f° r Na- 



aim to fit the raw experimental data with a minimal 
number of free parameters. We use the known crystallo- 
graphic structure in order to create the simulated crystal 
model[l5|. We then find all possible muon sites in the 
sample using electrostatic potential calculations In- 
deed the model finds 4 distinct muon sites, which was 



also seen previously in the pure NVO[13(. The next step 
was to find the correct spin configuration which results 
in the fields which give the muon precession measured in 
each site. Hence, dipolar fields were calculated with dif- 
ferent independent alignments of the AF moments of the 
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two zigzag chains, residing in the dc plane. The simula- 
tion also probed different IC helical modulations, that is 
the magnetic moment size at the I th unit cell is 



too cos(27rfc • I) 



(5) 



where k is the propagation vector along b and too is 
the V moment size, which is the only free parameter, 
and is given in Ref. [ill Q3 ■ We average over degener- 
ate states obtained by the simulation in order to accu- 
rately describe this model. We use \ 2 criteria as a crude 
guide in selecting the potential realized state. For com- 
parison, we also performed the same calculation for a 
linear IC spin density wave, which is proposed by neu- 
tron measurements [14J. However, as seen in Table I, x 2 
for a linear IC-SDW model was always higher than a 
helical-modulated setup. Figure EJa) displays the exper- 
imental ZF-spectrum of NVO taken at T = 2 K with the 
simulated spectrum assuming both the helical and linear 
configuration. 

We find that the best fit to the experimental data 
is, thus, a configuration with an helical IC modulation 
along b axis. Table Q] summarizes the simulation results 
and shows the propagation vector k versus x and the 
X 2 goodness-of-fit to the T = 2 K experimental data 
of these compounds. The value of NVO is found to be 
very consistent with that obtained by the recent neu- 
tron experiment [13], namely, the modulation period is 
5.214 x b along the 6-axis. However, the powder neutron 
diffraction result does not suggests a helical IC order but 
rather a linear IC-SDW order in NVO. The reason for 
this discrepancy is not clear. We speculate that the since 
the neutron result was obtained from powder samples 
through the Rietveld analyses, and the result published 
may be that of a local minima. Secondly, it might be 
that the dipolar field alone, which do not take into ac- 
count crystal fields might not be fully appropriate to de- 
scribe such a subtle magnetic state. In fact, a helical IC 
ordered state is more reasonable in explaining the coexis- 
tence of the AF order and metallic conductivity below Tm 
for NVO. Hence, /i + SR is found to play a significant role 
to determine the AF structure in NCVO. In Fig. EJb), 
we plot the corresponding magnetic structure suggested 
by the /i + SR simulation. 



B. Phase diagram 

To summarize the muSR measurements of NCVO car- 
ried out in the present work, we propose in Fig. [7] the 
complete phase diagram as a function of x. In partic- 
ular, the IC-AF phase with its wide field distribution, 
is present at low temperature of < x < 0.66. Interest- 
ingly, (above the IC phase), the compounds < x < 0.41 
show a mixture of CVO and the x — 0.416 compound by 
the TF measurements (shown as mixed- AF in Fig. [7]), 
but such mixture behavior is absent in the ZF data, in 
which only the IC phase is identified. The metallic phase, 
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FIG. 7: (Color online) The phase diagram of Na^Cai-a; V2O4 
family. PM: Paramagnetic, AF: Antiferromagnetic, IC: In- 
commensurate. 



x > 0.83, is characterized by numerous AF phases. AF1 
indicates the AF order identified by a single frequency, 
found below Tjy for x — 0.83 and at higher tempera- 
tures for x > 0.875. AF2 and AF3 describes the low 
temperature AF order characterized by 4 frequencies oc- 
curring in the x > 0.91. We therefore updated the previ- 
ously published phase diagram fl3| with the introduction 
of x = 0.41, identified the IC phase at x < 0.66, as well 
as the new data 0111 = 0.875, 0.958 shed additional light 
on the multiple AF phases. 



V. SUMMARY 

In summary, /iSR measurements on the 
Na a; Cai_ 2 ;V204 family reveal a diverse number of 
magnetic phases for different x. Compounds in the 
insulating phase, < x < 0.8, indicate a single muon 
site, data analyses identifies an incommensurate spin 
density wave phase. Whereas in the metallic phase, at 
high Na concentrations, x > 0.8, data suggests several 
muon sites. Numerical simulations verifies 4 magneti- 
cally distinct muon sites, and validates the existence of 
an incommensurate helical spin density of the V zigzag 
chain, with a spin-wave propagation vector k along b 
axis. Finally, we wish to point out that, although /x + SR 
is usually thought to be incapable to determine the 
correlation length of the magnetic order, the present 
study clearly demonstrates that /i + SR plays a crucial 
role for clarifying the magnetic structure of NCVO. 
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